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The relationship between the dimensionality of the cosmic ray pro- 
pagation function and the statistical distribution of anisotropies is 
demonstrated and an argument is presented in favor of "essentially" one 
dimensional propagation. This implies that a fluctuation explanation of 
low observed anisotropy can not be ruled out as has been stated by previous 

authors . xxxxxxx xx xxxxx xx xx xx xx xxxx x xxx xx x 

In a recent letter Ramaty et al 1 report results of a Monte Carlo cal- 
culation in which the injection of cosmic rays into the galaxy is considered 
to be a sequence of random discrete events in space time. In their letter 
they state that their results are in conflict with a suggestion 2 of the 
present author that small values of the cosmic ray anisotropy could result 
from the statistical nature of the injection mechanism. This remark is 
based on the result 3 that in the distribution of anisotropies small values 
are suppressed and the maximum liklihood value is of the order of the R.M.S. 
value. Ramaty et al also state 1 that a fluctuation origin of small anisotropy 
is only possible in the case of "strictly one dimensional" propagation of 
cosmic rays in the galaxy. This is demonstrated 3 by the result that when 
two of the dimensions were suppressed in the Monte Carlo calculation the 
suppression cf small values was not observed and a distribution was obtained 
that war. flat down to zero. 

It is the purpose of this note to point out the reason for this dependence 
on dimensionality, to demonstrate that small values of anisotropy are possible 





if the propagation is "essentially" one dimensional on a scale of the 
order of the Larmor radius of a cosmic ray particle, and to argue that 
this, in fact, is the case in our galaxy. In the following we shall con- 
sider the hulk cosmic ray flux density rather than the anisotropy & since 
the former has simpler statistical properties and the latter quantity is 
simply related to it by <jT where is the cosmic ray 

particle density. 

In the one dimensional case we note that the flux is the sum of many 
contributions from different injection events (supernovae explosions, pulsars 
etc.). Therefore, by the central limit theorem the distribution of T will be 
approximately Gaussian 

Pa ) dt ~ (t r<J *}) ' e*p(~J /?<- r '£) JT <n 

We shall assume a Gaussian form in the following argument since the results 
should be qualitatively the same for reasonable deviations from this form. 

This distribution is seen to be flat in the vicinity ofJ""= O . However, 
in the three dimensional case J is a vector with a distribution 


=( try'" {<!,'> <T*l<T,'>J 


and if < j A ' ^ *1 the quantity J */,£/ I s distributed as 
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The origin of the suppression nearj* *» O is now evident, it is the factor 


T* 'in the phase space volume element //y J 
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If on the other hand, instead of the R.M.S. values of the components 
being equal we have (0* 4B =/ )> straightforward integration 

gives a distribution for ^ j of 

P(T)dj ~ (2/ /r ) (fix*) exp(-5}%«) £ />// (T 8) J AT ( 4 ] 

where ( X ) — -i&rfCiX ) and £ - [ (**} J 

Using the asymptotic expression 4 , &rfi(x) —? A £Ap (X^) 

for XPO and — 7* <2Xp(x X )/&X as X we have 
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Clearly if <X = ft > B = O and the first case is always applicable; 

expression ( 5 ) becomes identical to expression ( 3 ). On the other hand if 

Bz(ltX) and the distribution is essentially equal to 

expression (l) (i.e. one dimensional) as long as J PP (%P.) 

We must now inquire as to what type of propagation of cosmic rays can 

bring about a distribution for the bulk flux of the form of expression (2) 

with the condition <( Tx ^ ^ <( ( . If the 

current ^ is the sum of many small currents from individual injection events 
/v * 

i.e. j ~ 22 with /V large it can be shown 5 that the dis- 


, tribution of £ approaches the form 


ihVin KNiinsHtHMttt 



- 4 - 


p(£)J l T - (irf 3 INI ^ ^ J >* ] J 3 J 


( 6 ) 


where J = jt& v an&(Al ) etn is the inverse 

of the covariance matrix M lm - NCjt j™ Xv , IN) I being 

its determinate „ In the above expressions^" ^ denote averages over the 
statistical variables 2 . Expression (6) is not yet in the form of (2) but one 
should note that J^f (and hence ) is a real, symmetric matrix and that 

it may always be diagonalized. If we can now find a propagation function 
that will give average values of ^ and that are small ( (9( e) J compared 
to j ^ we can have 

m. - <fii *<7</f i = ^ -- #?*> 


and ,/ j -O(l) then it is straightforward to verify that the values of Mu 
for the diagonalized matrix (the eigenvalues) will have the proper ordering, 
one of 1 ) and two of $C<: x ) . 

If we first consider three dimensional diffusion in a homogeneous magnetic 
field we may note while the mean free path along the field can be essentailly 
any value depending on the density of small irregularities in the field the 
diffusion mean free path in a direction perpendicular to the field is limited 
to be of the order of the Larmor radius P>L ■ “ we therefore choose 

"k (7^ and J we have upon performing the suitable 

averages 1 ? 2 (we here ignore the complication of the time dependence which 
while having problems of its own 2 adds essentially nothing to the ordering 
problem) 
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J 04 V<5 5 > etc. 
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with obvious permutation of X, ^ and ^ . 

This gives the right ratio of the various averages but absolutely it 
gives <y/^> = £(S ) «n a ( ~ &('/* * ) a very 

large ^ rather than a small value for and (fa ^ . In 

fact we have OCT £ O/**) 

where jt? denotes the cosmic ray density produced by an injection event and 
we note that this comes about because in this model the cosmic rays diffuse 
along a thin flux tube never getting more than a few times away from the 
line of force threading the injection event. If the cosmic rays are very 
lumpy (or stringy) we would expect this large fluctuation in the density. 
However, on the basis of meteorite data 6 we may rule out this model. 

On the other hand, if we consider the effect of the ergodlc nature of the 
galactic magnetic field 7 we see that cosmic rays will be transported across the 
mean field in a random or diffusive manner. The correct transport equation for 
this situation would describe diffusion in ^ along the mean field in time 
with diffusion in the X and directions with the position along the axis 
playing the role of time 7 . The diffusion coefficients are of the form 
/b ir and /& 2- . As a rough approximation to the more 

correct situation one can consider simply three dimensional diffusion with 
equal diffusion coefficients as was done by Ramaty et al. Hoover, it should 
be noted that, like all diffusion theories, it only applies on a large scale. 
The characteristic length in this case is the correlation length L of the 
random part of galactic field and L ~ to & f c In striking contrast 
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the distance scale that is applicable in any measurement of cosmic ray 
flux is the Larmor radius of a typical cosmic ray particle /? fc ^ 3^/0 p C . 
Therefore in computing local current densities we should consider the field 
to be smooth and homogeneous. If one treats the large scale transport of 
cosmic rays as an isotropic diffusion process with a mean free path /( the 
observed currents will have one component along the mean field with a R.M.S. 
value typical of /? but and ^ will be smaller by 

a factor Au//J ' 

In this situation our previous analysis applies and the observed flux 
will be distributed with a quasi Gausian with / assume 

for values of J" below this value. 


'A * '*»' 

for values of J { J ^ and wl ‘ th a rather sharp cutoff 


For any type of propagation function there will be a characteristic 

length/^ and a characteristic time . We have 2 by simple dimensional 

- . . 1 ^ _ 

analysis Al~h(i-c J where h is the injection event rate per unit volume 

"* y /i . j _ ^ v . .*■ . ** / 


and (f ^ ^ • This S^es { j ^ ~ fr L c % 

diffusion L c ~( j A C 7c) ’ ^ we assume with Ramaty et al that - 3 C 


and for 
* . 


then 


L c ~ 'Z.JToOfC and (J ^ 


It will be shown in a future publication that the streaming velocity is 
distributed with a singular distribution [P C'V) ^ W for large , 
hence there are no moments] that has a characteristic scale 0^. 

where (Sp^s is the mean square 

fluctuation of the cosmic ray density. Since 2 /V ^ 

We note that the median and J.O~ and ?CT levels of the anisotropy reported 
by Ramaty et al are proportional to /] for /I 7 p C. (large enough 
for nuclear collisions to be neglected). If we take their median value W 
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to be a good measure of (F^ our analysis indicates that the dis- 
tribution of £ should be quite flat down to values of about S //L ) 

Fz JO independent of /I or far smaller than anything observed. 

-y 

The point of all this is simply that although a value of & 4 
may be very unlikely with a flat distribution of this type it is no more unlikely 
than any other value measured with a precision of « lo’ 4 . From this 
we see that a small value is no more inconsistent with ^ ? / ft. than 
any other value. Furthermore it could be pointed out that if one chooses 


a constant ^ / O * years as do Jokipii and Parker 7 one obtains /) 




and a median £ of ~ 10*" 4 would not be obtained until pc 

In addition one could argue against simple diffusion from a point 
source by noting that the estimated 8 energy released in a typical injection 
event ~ 10 51 ergs is of the order of the ambient energy (a few eV per cc) 

n j 

contained in /0 (f>i) . This would indicate that each event would violently 
disrupt the galaxy at least over distances of the order of its thickness 
~ 200 pc. Furthermore if the non-linear propagation equation proposed by 
Skilling 9 is correct it would appear that the proper transport equation 
would more closely resemble the simple function originally employed 2 by the 
present author than a diffusion function. 

Perhaps the most serious argument against using the ergodic field line 
concept for Justifying the three dimensional transport of cosmic rays is the 
fact that it will not, in fact, rapidly disperse particles that are injected 
on neighboring field lines 7 . If a given field line is dispersed as 

z, ( m/g, x ) C tx (p) -'*<•« C XA (&) 

the correlation function for the random part of the magnetic field 

Z) A simple extension of the argument of Jokipii and Parker 7 

gives for two field lines separated by a distance at Jr « 0, 

((to A, /& t X (2L/8*) [Cjoi(o) ~ C xx (£)J , 
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For a Gaussian correlation function C xX (^)o< g.X />(-??*/} L % ) and we 

have ('(A*i ^ /t jj* 

or the relative dispersion of two neighboring lines is a factor (9//L ) smaller 
than the total random wandering of a given field line. Since /_ it /C& j 3 d. 
or of the order of the galactic disk thickness, a source of size /J? much 
smaller than this will produce a tangled thread of cosmic rays but not a 
diffuse cloud. 

All in all it would appear much to early to rule out the possibility 
that a low value for the observed cosmic ray anisotropy might be just an 
accident of our particular position in space and time. 
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